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ABSTRACT 

The  treatise  seeks  to  find  the  optimum  shape  of  a  thin 
beam  under  the  influence  of  lateral  buckling.   The  specific 
formulation  is  by  means  of  the  isoperimetric  problem  of  the 
calculus  of  variations.   An  energy  method  approach  yields 
the  governing  equations  for  the  problem.   An  analytic 
solution  was  not  obtained  due  to  the  nonlinearity  and 
coupling  conditions  of  the  equations. 

A  variable  height  lateral  buckling  problem  is  formulated 
as  an  alternate  attempt  to  find  the  optimum  design.   Through 
Bessel  equation  analysis  numerous  designs  are  obtained  and 
the  resulting  buckling  loads  are  calculated.   The  largest 
buckling  load  corresponds  to  a  beam  design  which  appears  to 
be  very  close  to  the  optimum  shape. 
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I.   INTRODUCTION 

A.   THE  OPTIMIZATION  PROBLEM 

The  performance  of  an  elastic  body  is  measured  by  the 
maximum  load  which  the  body  is  able  to  withstand.   The 
optimization  problem  considered  in  this  treatise  is  that 
of  determining  the  specific  design  of  an  elastic  structure 
such  that  the  structure  will  exhibit  a  maximum  performance 
capacity.   Herein,  the  term  design  refers  to  the  actual 
physical  distribution  of  material  within  an  elastic  body. 

The  optimization  problem  is  formulated  by  the  calculus 
of  variations.   This  method  offers  certain  advantages  in 
the  development  of  a  theory  for  optimization  of  an  elastic 
body.   The  specific  formulation  of  the  problem  is  by  means 
of  the  isoperimetric  problem  of  the  calculus  of  variations. 
By  this  method  the  problem  becomes  an  extremization  of  the 
total  potential  energy  of  the  system  with  a  superimposed 
volume  constraint.   The  extremization  of  the  functional, 
augmented  by  Lagrangian  multipliers,  leads  to  a  set  of 
Euler  equations  which  govern  the  design  of  the  structure. 

The  calculus  of  variations  is  no  new  technique  for 
obtaining  an  optimum  design.   Lagrange  [1]  formulated  an 
optimum  column  problem  by  this  means  in  1773.   He  sought  to 
determine  the  shape  of  a  column  for  which  the  buckling  load 
was  largest  among  all  columns  of  given  length  and  volume. 
Although  he  failed  in  obtaining  the  correct  solution,  his 


recognition  of  the  problem  belonging  to  the  calculus  of 
variations  inspired  others  to  continue  research  in  this 
field.   Clausen  [2]  later  provided  a  solution  to  this  prob- 
lem for  columns  of  circular  cross  sections.   Keller  [3] 
recently  showed  that  the  strongest  column  is  not  one  of  a 
circular  cross  section,  but  instead,  is  one  which  has  an 
equilateral  triangle  as  a  cross  section. 

The  isoperimetric  problem  for  optimum  design  received 
an  additional  stimulus  from  Keller  and  Tadjbakbsh  [4] .   This 
work  provided  a  variational  technique  for  minimum  volume, 
elastic  design  of  columns.   Taylor  [5]  treated  the  same 
problem  by  means  of  an  energy  method  approach.   In  his  work 
Taylor  showed  that  the  same  governing  equations  may  be 
obtained  by  direct  extremizati-on  of  a  functional,  which 
consists  of  the  system  potential  energy  and  a  volume  con- 
straint.  Taylor's  approach  is  the  essential  method  used 
for  optimization  in  this  thesis.   The  method  is  particularly 
useful  since  it  provides  a  simple  and  direct  means  for  iso- 
perimetric problem  formulation. 

An  extension  of  the  Taylor  approach  was  accomplished  by 
Salinas  [6]  and  provides  some  conditions  for  the  validity 
of  the  potential  energy  formulation.   Additionally,  the 
author  demonstrates  the  equivalence  of  maximum  load,  minimum 
volume  optimum  designs  and  presents  a  general  formulation  for 
optimization  of  elastic  structures. 
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B.   SCOPE  OF  THE  THESIS 

Using  the  Taylor  approach,  the  optimization  problem  is 
discussed  and  developed  by  use  of  Lagrangian  multipliers 
and  the  calculus  of  variations.   The  resulting  set  of  Euler 
equations  are  seen  to  be  the  system  equilibrium  equations 
and  an  optimality  condition.   Example  problems  are  included 
to  illustrate  the  usefulness  of  the  Taylor  method.   These 
include  a  strength  problem  for  a  uniformly  loaded  beam  and 
an  optimum  design  problem  for  a  particular  case  of  column 
buckling. 

The  method  is  then  extended  to  the  case  of  a  thin  beam 
subjected  to  lateral  buckling.   Under  lateral  buckling,  the 
beam  is  considered  to  buckle  in  a  plane  perpendicular  to 
the  vertical  plane  of  loading.-  The  buckling  configuration 
consists  of  rotation  and  sidewise  bending  away  from  the 
horizontal  axis  of  symmetry.   The  resulting  equations 
governing  an  optimum  design  are  found  to  be  non-linear.   In 
addition,  the  equations  are  coupled  in  rotation  and  bending 
All  attempts  for  an  analytic  solution  proved  futile. 

An  attempt  to  find  the  optimum  shape  has  been  made  by 
assuming  geometrical  conditions  for  instability  and  design 
of  the  beam.   A  variable  height  problem  has  been  formulated 
by  means  of  a  method  described  by  Federhofer  [7] .   Numerous 
designs  are  obtained  by  use  of  Eessel  equation  analysis  and 
a  predicted  optimum  shape  is  forseen. 

Throughout  the  discussion  all  problems  are  considered 
only  for  the  elastic  range  of  the  material.   Additionally, 
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the  material  is  considered  to  be  homogenous  and  isotropic 
All  problems  analyzed  are  restricted  to  rectangular  cross 
sections . 
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II.   PROBLEMS  OF  OPTIMAL  STRUCTURAL  DESIGN 

A.   THE  BASIC  PROBLEM  OF  STRUCTURAL  MECHANICS 

The  classical  theory  of  structural  mechanics  is  concerned 

with  the  behavior  of  continuous  elastic  bodies  which  are 
subjected  to  specific  loading  conditions.   In  general  the 
components  which  define  mechanical  behavior  are  the  displace- 
ment vector  u,  and  the  stress  and  strain  tensors,  a.,  and 

The  strain  tensor  is  explicitly  defined  by  the  geometry 
of  deformation  of  an  elastic  body.   Hence,  the  strain- 
displacement  relation  may  be  represented  by 

,   §_  =  e(u)  -  (2.1) 

From  the  constitutive  laws  of  elasticity  it  is  known  that 
stress  is  related  to  strain,  that  is 

a  =  a(e)  (2.2) 

Thus,  stress  and  strain  are  both  representable  as  functions 
of  the  displacement.   In  classical  theory  any  one  of  the 
three  quantities,  u,  a,  and  e  is  sufficient  for  determining 
the  other  two. 

The  mechanical  state  of  a  structure  is  defined  by  a 
boundary  value  problem  which  is  comprised  of  equilibrium 
equations  and  the  associated  boundary  conditions.   Thus,  a 
classical  problem  in  structural  mechanics  consists  of  using 
the  constitutive  relations  between  stress,  strain  and 
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displacement  to  derive  and  solve  the  governing  equilibrium 
equations  and  their  associated  boundary  conditions. 

B.   CALCULUS  OF  VARIATIONS 

Problems  concerned  with  the  determination   of  extreme 
values  of  integrals  whose  integrands  contain  unknown  func- 
tions below  to  that  field  known  as  the  "calculus  of 
variations"  [8] .   Such  an  extremal  refers  to  the  maximum, 
minimum  or  stationary  value  of  a  functional.   A  functional 
is  a  real  valued  function  of  a  function. 

Only  the  integral  functional 

I[f(x)]  =  /  F[f(x)]dx  (2.3) 

D 

will  be  considered  herein,  where  D  is  the  n-dimensional 
domain  over  which  the  n-dimensional  vector  x  ranges.   The 
function  which  extremizes  I[f(x)]  is  denoted  as  f*(x)  and 
the  associated  value  of  the  functional  is  I*  (i.e.  I*  = 
I[f*(x)]). 

In  general,  problems  in  structural  mechanics  concerning 
optimization  will  assume  the  form 

/  F(x;u,u'  ,u", <»,<}>'  ,.  .  .)dx  (2.4) 

Rl 

for  a  one  dimensional  case.   A  one  dimensional  case  is 

generally  characterized  as  having  one  dimension,  its  length, 

much  larger  than  the  other  dimensions.   In  the  above  form  x 

represents  an  independent  variable.   The  quantities  u,  (J)  and 

their  respective  derivatives  represent  dependent  variables. 

In  this  thesis  u  and  <f>  represent  displacement  functions  or 
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state  variables.   The  nature  of  the  specific  problem  governs 
the  number  of  dependent  and  independent  variables. 

Performing  a  variation  with  respect  to  a  functional  is 
denoted  by  the  symbol  6.   From  [8]  it  can  be  shown  that  a 
stationary  value  of  the  functional  occurs  for  6F  =  0.   This 
operation  yields  what  is  known  as  a  set  of  Euler  equations. 
For  example,  if  equation  (2.4)  is  represented  by  the  par- 
ticular form 


I  =  /  F(x;u,u'  ,u",<t>,cj>f  ,A)dx 
D 


(2.5) 


then  performing  the  extremal  operation  yields  three  Euler 
equations  [91 • 


5    F   =    0 :      r—  -   tc— 

U  dU  dX 


6    F   =    0:       H  -  -r- 

<J>  d<J>  dX 


AF   "    0:      fl=    ° 


f3F 

+ 

9 

8F    | 

l3u'. 

9x2 

13u" 

i          / 

fBF 

= 

0 

=    0 


(2.6) 
(2.7) 
(2.8) 


plus  a  set  of  natural  boundary  conditions. 


C.   THE  ISOPERIMETRIC  PROBLEM 

The  isoperimetric  problem  in  the  calculus  of  variations 
provides  for  extremizing  functionals  whose  dependent  vari- 
ables are  constrained  by  an  integral  relation.   The  integral 
constraint  is  accommodated  by  constant  Lagrangian 
multipliers . 

As  an  example,  suppose  that  the  maximum  (or  minimum) 
value  of  an  integral 
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I  =  /  F(x;u3u' )dx  02.9) 

R 

is  sought,  subject  to  the  condition  that  another  integral 

J  -  /  G(x;u,uf )dx  (2.10) 

R 

is  known  to  have  a  constant  value.   Using  the  method  of 

Lagrange  multipliers  [8]  an  augmented  functional 

I+XJ  =  /  [F(x;u,uT)  +   G(x;u,u')]dx  (2.11) 

R 

may  be  constructed  and  its  free  extremum  considered.   The 
necessary  condition  for  a  stationary  value  is  the  Euler 
equation 

3(F+  G)    _8_   9(F+  G)  ,  . 

8u      8x     3u'  \*.±*) 

where  F  has  been  replaced  by  F-*  =  (F+AG). 

D.   POTENTIAL  ENERGY 

The  total  potential  energy  of  a  conservative  elastic 
system  can  be  expressed  as 

T  =  U  -  V    .  (2.13) 

U  represents  the  total  strain  energy  of  the  structure,  and 
V  is  the  potential  energy  of  external  forces  associated  with 
the  displacement  of  the  system  under  any  set  of  given  load- 
ing conditions.   The  total  strain  energy  of  the  system  is 
the  volume  integral 

U  =  \  I  o,,e. ,dV  (2.14) 

2  v   ij  ij 

where  a.,  and  e        are    components  of  the  stress  and  strain 
U       ij 

tensors . 
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The  potential  energy  of  external  forces  is  defined  as 
the  inner  product  of  the  external  force  vector  P  and  the 
associated  surface  displacement  vector  u,  that  is, 

V  =  /  P  •  u  ds 
s 

where  s  is  the  domain  over  which  P  is  defined. 

E.   PRINCIPLE  OF  MINIMUM  POTENTIAL  ENERGY 

A  powerful  tool  in  the  field  of  structural  mechanics  is 
the  principle  of  minimum  potential  energy.   This  theorem 
states  that  the  total  potential  energy  of  a  conservative 
elastic  system  attains  a  minimum,  with  respect  to  variation 
of  kinematically  admissible  displacement  functions,  when 
the  displacement  function  satisfies  equilibrium. 

The  total  potential  energy  has  previously  been  defined 
in  terms  of  strain  energy  and  external  work.   These  quanti- 
ties may  be  expressed  as  functions  of  the  displacement 
vector  u.   Thus,  the  equilibrium  equations  and  associated 
boundary  conditions  are  derivable  from  a  variation  with 
respect  to  the  displacements  of  a  potential  energy  functional 

The  variational  operation 

6U(U-V)  =  0  (2.16) 

can  then  be  performed  for  which  the  resulting  Euler  equations 
physically  represent  the  system  equilibrium  equations.   The 
optimization  method  considered  in  the  thesis  is  based  on 
this  important  principle. 
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F.   A  STRUCTURAL  OPTIMIZATION  PROBLEM 

As  an  example  of  structural  optimization,  consider  the 
problem  of  obtaining  the  optimum  shape  for  column  buckling. 
The  problem  may  be  characterized  by  the  statement :   "For  a 
column  of  given  length  and  volume  of  material,  determine  the 
column  shape  for  which  the  Euler  buckling  load  is  a  maximum." 

The  column  shape  is  considered  to  be  such  that  all  cross 
sections  are  similar.   The  moment  of  inertia  I(x)  is 
assumed  to  be  related  to  area,  A(x) ,  that  is 


n 


EI(x)  =  kA  (x) 


(2.17) 


where  E  represents  Young's  modulus  and  k  is  a  constant 
depending  on  the  section  considered.   For  example,  the 
moments  of  inertia  of  rectangular  sections  of  variable- 
width-constant  height,  and  variable  height-constant  width 
are  (Fig.  1) , 


Figure  1.   Column  Shapes  Analyzed,  (a)  variable  width  - 
constant  height,  (b)  variable  height  - 
constant  width. 
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EIyy  =  Eh  ^X)  =  kA(x)  (n  =  1)        (2.18) 

for  case  (a)  and 

EIyy  =  Ebh^x)  =  kA3(x)  (n  =  3)        (2.19) 

for  case  (b ) . 

For  the  Euler  buckling  problem  the  total  potential 

energy  is 

I 
T  -  \  S       (EIv"2  -  Nv'2)dx  (2.20) 

d    0 

where  v  represents  lateral  displacement  for  the  buckled 
state.   Introducing  an  isoperimetric  condition  for  a  con- 
stant volume  of  material  V  gives 

o 

/  A(x)dx  =  V  (2.21) 

0  ° 

The  volume  constraint  is  accommodated  through  the 
formation  of  the  functional  T* ,  defined  by 

T*  =  T  -  AV  (2.22) 

o 

where  A  is  a  constant  Lagrangian  multiplier.   Combining 

(2.17),  (2.20)  and  (2.22)  the  constrained  functional  is 

I 
T*  =  /   (-  An  v"2  -  ^  v'2  -  AA)dx  .  (2.23) 

0    2  2 

The  governing  equations,  in  addition  to  the  isoperimetric 
equation,  are  obtained  by  performing  a  variation  with 
respect  to  v  and  A  to  obtain 
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6yT*  =  0:   (kAnv")"  +  Nv"  =  0  (2.24) 

6AT*  =  0:   |nkAn_1  v"2  "  X  "  °  (2.25) 

Equation  (2.24)  is  the  equilibrium  equation  and  (2.25)  an 
optimality  condition. 

The  non-linear  equations  (2.24),  (2.25)  and  (2.21)  define 
the  optimum  design  with  respect  to  constant  volume  -  maximum 
load.   In  general  these  are  extremely  difficult  to  solve 
in  closed  form.   However,  an  exact  solution  can  be  obtained 
for  case  (a)  where  n  =  1.   These  governing  equations  then 
become 

(kAv")"  -  Nv"  =0  (2.26) 

|  kv"2  =  X  (2.27) 

£ 
/   A(x)dx  =  V  (2.21) 

0 

Equation  (2.27)  implies  constant  bending  curvature,  i.e. 

v"  =  /2X/k  (2.28) 

Substituting  (2.28)  into  (2.26)  and  integrating  twice  gives 

M      2 

A  =  P-   (2_  +  cnx  +  CQ)  (2.29) 

k   2     1      2 

Similarly 

2 

v  =  /2X/k   (~-  +  C3x  +  C^)  .  (2.30) 

If  the  column  is  considered  to  be  simply  supported,  the 
boundary  conditions  for  displacement  are 

v(o)  =  v(£)  =  0  (2.3D 
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£ 
for  which   C^   =   -  r,    Cn    =   0,    and 


x 


v(x)    =    /2X/k       (^-  -   ^f)     .  (2.32) 

The   moment    conditions    are 

M(o)    =   M(£)    =    0  (2.33) 

I 
from  which  C,  =  -  r,  C      =  0,  and 

2 
A(x)  =  |  (^-  -  ^)  .  (2.34) 

The  critical  load  is  now  determined  from  the  isoperimetric 
condition 


for  which 


_12kV0    -Eh2V0 


N   .  .     =  5-^  =  =-^  .  (2.36) 

optimum     ^3       £3 

Substituting  this  expression  into  (2.3*0  determines  the 
optimum  design  in  terms  of  known  constants,  that  is 

V  0 

A(x)  =  — -2-  (ix   -   xd)    .  (2.37) 

121 J 

Note  that  the  Lagrange  multiplier  does  not  appear  in  either 
of  the  expressions  for  A(x)  or  N.   This  is  in  distinction 
to  strength  optimization  problems  where  X  is  a  representa- 
tion of  maximum  strength  and  remains  in  the  final  analysis 
(see  Appendix  A). 
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III.   THE  LATERAL  BUCKLING  PROBLEM 

A.   GENERAL  DESCRIPTION 

The  problem  of  lateral  buckling  is  important  in  the 
design  of  beams  which  have  no  lateral  support  and  possess  a 
shape  such  that  the  beam  height  is  large  in  comparison  to 
the  width.   Such  a  configuration  may  become  laterally 
unstable  when  the  loading  reaches  a  critical  value.   The 
resulting  failure  or  "buckle"  is  a  combination  of  lateral 
bending  and  twisting  away  from  the  horizontal  axis  (Fig.  2) 
Prior  to  buckling  the  load  needed  to  cause  this  failure  may 
be  considerably  less  than  the  corresponding  load  that  would 


Figure  2.   Simply  Supported  Beam  In  a  Lateral  Buckled  State 
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cause  bending  failure.   This  is  due  to  a  large  flexural 
rigidity  of  the  beam  in  the  vertical  plane  of  bending  as 
compared  to  a  much  smaller  lateral  bending  resistance  in 
the  plane  perpendicular  to  the  applied  load. 

A  brief  but  informative  discussion  of  the  historical 
development  of  the  lateral  buckling  problem  is  given  by 
Bleich  [12]  and  will  not  be  included  here.   The  methods 
used  for  development  of  the  problem  in  this  thesis  are  sim- 
ilar to  those  used  by  Sechler  [13]  and  Timoshenko  [14] . 

B.   DEVELOPMENT  OF  CLASSICAL  SOLUTION 

While  the  solution  to  the  classical  lateral  buckling 
problem  may  be  found  in  a  number  of  texts,  its  inclusion 
herein  is  considered  necessary  for  understanding  the 
descriptive  geometry  of  the  problem  and  later  development 
of  the  problem  of  a  variable  height  solution.   The  example 
to  be  considered  is  that  of  a  simply  supported  beam  loaded 
at  its  centroid  by  a  concentrated  load  P  at  the  origin 
(Pig.  3).   The  centroidal  loading  condition  is  necessary  to 


Figure  3-   Simply  Supported  Beam  Loaded  at  the  Centroid 
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simplify  the  problem.   Loading  at  any  point  other  than  the 
centroid  introduces  an  additional  work  term  which  leads  to 
some  mathematical  difficulties.   This  condition  is  described 
in  Appendix  B. 

It  is  assumed  that  the  load  P  remains  in  the  vertical 
xy  plane  throughout  buckling.   The  ends  of  the  beam  are  so 
constrained  as  to  prevent  rotation  in  the  yz  plane.   The 
buckled  configuration  is  shown  in  Figure  4. 

The  deflection  of  the  beam  is  defined  at  any  point  x 
from  the  origin  by  the  components  u,  v  and  <J> .   Components 
u  and  v  represent  the  displacement  of  the  centroid  in  the 
z  and  y  directions  while  <\>    is  the  angle  of  rotation  of  the 
cross  section  in  the  yz  plane.   Thus,  the  boundary  condi- 
tions , for  the  problem  are 

uU/2)  =  4>U/2)  =  0  (3.1) 

-  u'(0)  =  <f>'(0)  =  0  (3.2) 

The  fixed  coordinate  axes  are  labeled  x,  y  and  z.   In 
the  buckled  state  the  principle  axes  of  symmetry  at  any 
section  mn  are  defined  as  £  ,  n ,  and  £.   The  usual  right 
hand  rule  convention  for  vectors  is  used  to  describe  posi- 
tive moments  with  respect  to  the  described  axes. 

By  the  usual  assumptions  of  small  deflection  theory  the 

equations  for  principle  moments  can  be  written 

H2n 
ET   £-§.  -  EI  V"  =  M  (3.3a) 


2k 
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Figure  k.      Lateral  Buckling  Configuration 
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H2r 

EI   2_k  ,  EI  u"  =  -M  (3.3b) 

n  d£.2    n      n 

GJ  ||  =  GJ$'  =  M?  (3.3c) 

These  moments  are  physically  described  as  follows: 

i)   M   is  the  usual  bending  moment  due  to  the 
5  load  P 

ii)   M   is  the  bending  moment  associated  with 
n  lateral  buckling 

iii)   M~  is  the  twisting  moment  associated  with 
lateral  buckling. 

With  respect  to  the  fixed  x,  y  and  z  axes  the  above 

moments  can  be  represented  [131  as 

Mr    =   M      +   M   u'    =    GJ<(>  '  (3.^a) 

C,  X  Ct 

M      =   M      +   M   4>      =    -EI    u"  (3.4b) 

n         y         z  n 

M      =    _m    <}>    +    M      =    EI  v"  (3.4c) 

c  y  z  ? 

Let  Cn  =  GJ  and  C0  =  EI  .   Differentiating,  the  first 
1  2     r|  ' 

equation  with  respect  to  x  and  substituting  for  u"  from  the 
second  results  in  a  general  equation  which  can  be  written 

M  2       ,  MM 

For  the  problem  considered  let  the  deflection  at  the 
beam  center  (maximum  deflection)  in  the  z  direction  be 
denoted  by  u, .   The  moments  at  any  distance  x  are  determined 


as 


Mx  =  -P/2  (u1-u)  (3.6a) 

M  =  0  (3.6b) 
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1VL  =  P/2  U/2-x)  (3.6c) 

Differentiating  and  substituting  into  equation  (3-5) 

yields 

2  £      2 
P  (p  -  x) 

4>"  +  "Tic^c *  =  °  '  (3*7) 

Introducing  the  new  variable  6  =  (£/2-x)  and  letting 

2 
2     P 
k   =  |,p  =    equation  (3.7)  becomes 

<}>"  +  k262<t>  =  0  .  (3.8) 

Equation  (3-8)  is  now  in  the  form  of  a  Bessel  equation  for 
which  the  general  solution  has  the  form 

<f>  =  Bh    [a^cf  32)  +  a2J^  (|  62)]  (3.9) 

where  J,  and  J  a  are  Bessel  functions  of  the  first  kind  of 
orders  \   and  -\   respectively  [1H] .   For  the  simply  supported 
beam  the  boundary  conditions  with  respect  to  c}>  are 

4>  =  0   at   6  =  0     (i.e.  x  =  1/2)  (3.10a) 

<t>'  =  0   at   6  =  1/2         (i.e.  x  =  0)  .  (3.10b) 

The  first  boundary  condition  yields  a?  =  0 .   Differentiation 
of  <$>   with  respect  to  6  now  yields 

||  =  alk63/2  J_3/1((|  B2)  .  (3.11) 

Applying  the  second  boundary  condition  shows  that  buckling 
occurs  when 

J_3/4  (-^~)  =  0  •  (3.12) 
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From  a  table  of  Bessel  functions  [151  the  first  root  of 
(3.12)  is 

k£2 

£g-  =  1.058508  .  (3.13) 

Substituting  for  K,  the  buckling  load  is  found  to  be 


16.936  V^i^I" 


Pcr  =  ,2  ^3.14) 

Equation  (3.1*0  represents  the  exact  solution  for  the  prob- 
lem considered  and  will  be  used  as  a  comparative  value  for 
critical  loads  calculated  for  variable  cross  sectional 
designs  later. 

C.   OPTIMIZATION  TECHNIQUE  APPLIED  TO  LATERAL  BUCKLING 

The  energy  formulation  described  in  Section  II  is  now 
used  to  develop  the  governing  equations  for  an  optimum  shape 
of  a  thin  beam  which  undergoes  lateral  buckling.   The 
assumptions  made  during  the  preceeding  section  apply.   Addi- 
tional assumptions  include  neglecting  the  strain  energy  due 
to  "warping"  [14]  and  elongation  of  the  beam. 

The  case  considered  is  that  of  a  variable  height-constant 
width  beam  (Fig.  5).   The  principle  moments  of  inertia  with 
respect  to  the  buckled  state  are 

b2 
In  =  ±  b3Q   h(x)  =  j|  A(x)  (3.15) 

b2 
J  =  |  b^  h(x)  =  -2-  A(x)   .  (3.16) 

The  strain  energy  due  to  lateral  bending  may  be  written 
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Figure  5.   Variable  Height,  Constant  Width  Beam, 


1/2 


„2 


U,  =  /     C,Au"   dx 
1    0      L 


(3.17) 


where  C,  =  Eb  /12  .   The  strain  energy  associated  with 
1     o  aj 

twisting  of  the  structure  is 


1/2 


,2 


Up  =  /     CpA<j)'   dx 
0 


(3.18) 


where  Cp  =  Gb  /3. 

As  the  beam  buckles  the  load  P  is  lowered  an  amount  which 
may  be  represented  as  [14] 


4>u"  {1/2   -   x)  dx 


(3.19) 


The  potential  energy  of  external  forces  associated  with  the 
buckling  configuration  is  determined  by 


1/2 
V  =  P  /     4>u"  {1/2   -  x)  dx  . 
0 


(3.20) 


The  total  potential  energy  of  the  system  is  U  -  V  or 

1/2 
T  =  /     [C  Au"2  +  CpA(f)'2  -  P4>u"  {1/2    -  x)]dx.  (3.21) 

To  formulate  a  functional  by  the  method  described  earlier, 
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an  isoperimetric  constraint  is  introduced  for  a  given  volume 

of  material, 

1/2 
f  Adx  =  V  (3.22) 

0  ° 

The  augmented  functional  may  be  written 

T*  =  /     [C  Au"2+  C  Ac})'2-  P<J>u"(4-x)  -  XA]dx.(3.23) 
0 

Performing  a  variation  with  respect  to  the  variables  u,  <t> 

and  A  yields  three  Euler  equations : 

£ 
<5uT*    =    0:       [2C1Au"    -    P<{>(§  -    x)]       =0  (3.24) 

6    T*    =    0:       202^')'    +    Pu"(|  -    x)    =    0  (3.25) 

6AT*    =    0:       C,uf|2    +    C^'2    -    X    =    0    .  (3.26) 

Note  that  the  variation  with  respect  to  the  state  variables 
u  and  <j)  yields  equilibrium  equations.   Variation  with 
respect  to  the  design  variable  A  introduces  an  optimality 
equation.   The  associated  boundary  conditions  are: 

u(|)  =  <D(|)  =  0  (3.27) 

u'(0)  =  <j)'(0)  =  0  (3.28) 

Integration  of  equation  (3.24)  twice  and  evaluating  at  the 
boundary  conditions  determines 

2C1Au"  -  P<J>(|  -  x)  =  0  (3.29) 

Integration   of    (3-25)    yields    the    result 

C2A<J>'    +   |  u'    (|  -   x)    +   |   (u1-u)    =    0  (3.30) 
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where  u,  represents  the  maximum  deflection  at  the  beam 
midspan. 

The  above  equilibrium  and  optimality  equations  along 
with  the  isoperimetric  condition  (3.22)  define  the  optimum 
shape  for  the  structure.   These  governing  equations  are  noted 
to  be  non-linear  and  an  attempt  to  solve  them  has  been 
unsuccessful.   Unfortunately,  these  equations  resist  uncoupl- 
ing, which  precludes  getting  to  the  isoperimetric  condition 
as  was  accomplished  easily  in  Section  III.   Additionally,  no 
computer  techniques  were  found  which  could  solve  systems  of 
simultaneous  nonlinear  equations. 

After  considerable  efforts  to  solve  the  above  non-linear 
system  proved  fruitless,  an  attempt  to  formulate  a  variable 
height  solution  from  a  specified  polynomial  was  undertaken 
in  hopes  of  predicting  or  converging  to  an  optimum  solution. 
The  method  used  and  results  obtained  are  included  in  the 
following  section. 
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IV.   A  VARIABLE  HEIGHT  BUCKLING  PROBLEM 


A.   DESCRIPTION  OF  THE  PROBLEM 

The  problem  of  lateral  buckling  of  a  variable  height, 
constant  width  beam  is  now  considered.   Again,  the  problem 
to  be  investigated  is  that  of  a  simply  supported  beam, 
loaded  at  the  midspan  centroid.   However,  the  height  now 
varies  according  to  the  power  law 


h  =  hQ(l  -  f) 


(4.1) 


where  h   is  the  height  at  the  center  section  of  the  beam 
o 

(Fig.  6).   The  value  of  n  is  restricted  to  lie  between  the 
interval  0  _<  n  <_  1  for  reasons  which  become  physically 
obvious  in  later  analysis  of  the  problem.   The  idea  for 
using  this  type-  of  formulation  was  based  on  the  results 


h=h0(l-2X/2) 


n 


Figure  6.   Variable  Height  Beam  Based  on  the 

Power  Law   h  (l-2x/£)n. 
o 
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obtained  for  a  variable  height  cantilever  problem  by 
Federhofer  [7] . 

Differentiating  equation  (3.4a)  with  respect  to  x  and 
noting  that  the  previous  constants  C,  and  C?  are  now  vari- 
able functions  of  x,  reduces  (3.4a)  to  the  form 

C,<|>"  +  C,  T<j>?  =  MT  +  M'u'  +  M  u"  (4.2) 

1T      1  Y      x     z       z 

Substituting  (3.4b)  into  (4.2)  for  u"  and  recalling  the 
moment  conditions  from  (3.6)  (and  their  respective  deriv- 
atives) reduces  (4.2)  to 

P2(*  -  x)2 
<T  +  ~  Cj_<t>'  +   h2c   c ♦  -  0  .  (4.3) 


Letting 


y  =  (1  -  Y}  (4'4) 


and  using  the  relations 

Cl  =  T2~  Et)3  h(x) 
C2  =  |  Gb3  h(x) 

the  lateral  and  flexural  stiffness  at  the  beam  center  may 
be  designated  as  C,  and  C«.   Thus,  the  new  variable  stiff- 
nesses can  be  represented  by 

C1  =  C°  (1  -  ^-)n  =  C°yn  (4.5a) 

C2  -  C°  (1  -  ^)n  =  C°yn  .  (4.5b) 

Using  the  chain  rule 
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d<f>    _   d<£      dx  fh    ^v 

d7  "  dl    d7  (4*6) 


and  noting 


x   =   |   (1-y) 

dx^   I*      dY 
Similarly 

_1_   iH  d^L  =  4__  n  d£.  fk   R} 

C1   dx   dx  '  ^2   y   dy   .  ^-o> 

Substitution  of  (4.7)  and  (4.8)  into  equation  (4.3)  yields 

^2,       ..    p204  2(l-n) 

d_|+n  di+^_j =Q_  (i|>q) 

dY2    Y   dY      64C°C° 

2" 
The  dimensionless  parameter  K  may  be  defined 

P2.4 

=  IT  .  (4.10) 


64C°C° 


Equation  (4.9)  may  now  be  written 

d£|+  n   dj>    K2  2(l-n)   =  (4  n) 

dY     Y    Y 

Transformation  of  equation  (4.11)  into  a  Bessel  equation  can 
be  accomplished  by  the  following  substitutions. 

P  =  ^  y2~n  (4.12a) 

1-n 
4>  =  Y   2   6(y)  (4.12b) 

Equation  (4.11)  is  then  reduced  to  the  form 
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dp    p   dp        p 
where  v  is  representative  of 

v  =  5T|fay  .  (lt.13a) 

Verification  of  the  transformation  is  accomplished  in 
Appendix  C. 

1 .   Non-integer  Value  of  v 

For  the  non-integer  values  of  the  parameter  v 
between  zero  and  one  equation  (4.13)  has  the  following  solu- 
tion [16] 

6  =  a1Jv(p)  +  «2J_V(P)  (4.14) 

where  J   and  J    are  Bessel  functions  of  the  first  kind  and 
v      -v 

order  v. 

Recalling  (4.12b),  (i.e.  $  -   y  '        3(y))3  the  complete 

solution  is  determined  by 
1-n 
4>  =  Y  2   [a1Jv(p)  +  «2J_v(p)]  (4.15) 

The  boundary  conditions  associated  with  the  lateral  buckling 
problem  are : 

<Ky  =  0)  =  0      (i.e.  x  =  |)  (4.16a) 

||  (Y  -  1)  -  0    (i.e.  x  =  0)  .  (4.16b) 

From  the  first  condition  and  expansion  of  J   ,  it  is  con- 
cluded that  cu  =  0 .   Therefore  (4.15)  becomes 
1-n 
*  =  Y  2   [o^J^p)]  (4.17) 
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for  which  differentiation  yields 


dy     1 


-1-n 


1-n 


1-n 


J  (p)  +  y    J'(p)  p- 
v  '      v  K   dy 


It  is  additionally  known  that 


dp_ 
dy 


=   Ky 


1-n 


and    [17] 


J'(P)  =         -        J  ...      (P)  +         -       J     ,,  N 

v    K  v+1  p      v(p ) 


(4.18) 


(4.19) 


(4.20) 


Evaluation  of  -r—   at  the  second  boundary  condition  with 
y  =  1  and  p  =  K/(2-n)  yields  the  solution 


j   (  K  )    u   ^K_ 


(Jkr) 


v  v2-n'    1-n   v+1  v2-n 


(4.21) 


The  buckling  loads  may  now  be  determined  by  solution  of 
equation  (4.21).  This  may  be  accomplished  by  independently 
calculating  each  case 


j   (  K  )  =  o 
v   2-n 


(4.21a) 


and 


K 


J 


(A)  =  o 


(4.21b) 


1-n   v+1  v2-n 

and  determining  the  initial  intersection  of  the  resulting 
Bessel  equation  curves.   For  each  considered  value  of  n,  K 
is  found  numerically  from  the  argument  of  J.   A  computer 


Handbook  of  Mathematical  Functions,  National  Bureau  of 
Standards  Applied  Mathematics,  Series  55  }  p.  363,  June 
1964. 
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solution  for  this  problem  is  included  in  the  Computer  Pro- 
gram section  of  the  thesis.   The  Bessel  equation  curves 
which  show  the  intersecting  roots  are  included  in  Appendix  D 
When  n=0  the  beam  design  reduces  to  that  of  the  rectang- 
ular case.   The  computed  solution  from  equation  (4.21a)  and 
(4.21b)  provides  an  argument  value  of  1.05851.   Thus,  the 
critical  load  is  determined  by 

=  1.05851  (4.22) 


2-n 


or 


16.936^0^0 


0^,0 


Pcr  =  J    -  2  (K. 23) 

precisely  the  value  determined  previously.   The  values  for 
other  solutions  are  discussed  in  a  later  section. 
2 .   Integer  Value  of  x> 

When  the  value  of  n  is  equal  to  one  the  resulting 
value  of  v  becomes  zero  and  equation  (4.13)  is  reduced  to 
the  form 

2 

dp2   p  dp 

The  complete  solution  of  this  equation  can  be  shown  to  be 
of  the  form  [17] 

6  =  a1JQ(p)  +  ct2yo(p)  (4.25) 

where  y  is  known  as  a  Weber  function  or  Bessel  function  of 
'  o 

the  second  kind.   Substituting  for  <J> ,  the  above  equation  is 
transformed  to 
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1-n 
<t>   =  y   2   [a1JQ(p)  +  a2yQ(p)]  .  (4.26) 

This  solution  is  valid  only  if  v  is  of  an  integer  order. 
Thus,  in  the  restricted  range  n  is  allowed  only  to  be  equal 
to  one,  and  the  y  multiplier  must  be  to  the  zero  power. 
Equation  (4.24)  reduces  to 

4>  =  a1Jo(KY)  +  a2yo(Ky)  .  (4.27) 

I 
The  value  of  y  (Ky)  for  y   =  0  (i.e.  x  =  j)    is  -°°  which 

implies  that  a   is  equal  to  zero.   Differentiating  cf>  with 

respect  to  the  remaining  term  and  evaluating  at  the  second 

boundary  condition  (y  =  1,  x  =  0)  requires  that 

-J1(K)  =  0  (4.28) 

for  the  solution  to  (4.26).   This  value  is  found  [181  to  be 
either  the  trivial  case  (i.e.  zero)  or  3-8317.   The  latter 
result  implies  a  singularity  condition  with  respect  to  the 
previous  non-integer  solutions.   A  graphical  representation 
(Fig.  7)  of  values  of  n  versus  the  roots  of  equation  (4.21) 
does,  in  fact,  imply  the  resulting  value  converges  to  zero 
as  n  approaches  one  as  a  limit.   Therefore,  the  trivial  case 
is  the  solution.   The  beam  is  apparently  unstable  for  any 
load  when  in  this  configuration. 

E.   CALCULATION  OF  AMENDED  LOADS 

Extension  of  the  problem  is  now  made  in  an  attempt  to 
find  the  design  among  all  those  analyzed  that  maximizes  the 
buckling  load.   Initially  h  was  described  as  a  fixed  height 
at  the  beam  midspan.   This  center  height  is  now  allowed  to 
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Figure  7.   Convergence  of  the  Foot  x  as  N  approaches  One. 

vary  such  that  each  variable  design  studied  will  have  the 
same  volume  as  the  constant  rectangular  beam. 

This  may  be  accomplished  by  assuming  a  constant  width 
for  all  designs,  integrating  the  power  law  (4.1)  and  setting 
the  resulting  value  equal  to  that  of  the  classical  rectang- 
ular beam.   Using  symmetry  and  integrating  the  power  law 

I 
from  x=0  to  x=  „- 


1/2  2x  n 

b  /  h   (1  -  =r)      dx  =  V 

o    °     £        ° 


(4.29) 


determines  the  volume  of  the  variable  design.   For  the  rec- 
tangular beam  the  volume  for  one-half  the  structure  is 

hrbd 

— p—  .   Thus,  equating  the  two  determines 

hQ  =  (n+1)  hr  (4.30) 

It  is  obvious  that  the  center  height  must  be  expanded  for 

the  variable  design  in  order  that  the  volumes  are  compatible 

When  h   is  restricted  to  a  fixed  value  the  critical 
o 

loads  may  be  calculated  as  described  in  Section  IV-A,  i.e. 
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using  the  intersecting  value  of  (4.21a)  and  (4.21b). 

p#,2 
Recalling  that  K  =  -—^   ,  the  critical  load  is  defined  as 

8  7^1 

O00 


8(2-n)xNyc°C 


= — x  2  (4  ?1) 

crit         %2  ^,Jlj 

where  x  represents  the  intersecting  value  of  the  Bessel 
equation  curves.  The  amended  load  is  now  found  by  sub- 
stituting for  C, ,  C0  and  h  . 

12      o 

8(2-n)xh 
Pcrit  =     £2    °  7b6  EG/36   ,  C4.32) 

Thus, 

8(2-n)x(n+l)h   

P*  ■  ~2 -J^T  ^-33) 

where 

C4  =  EG.  (4.35) 

The  amended  load  is  referred  to  as  P*  since  its  value  is 
greater  than  its  unexpanded  counterpart.   It  is  left  in 
the  above  form  in  order  that  it  may  be  readily  compared  to 
the  solution  for  the  rectangular  beam. 

A  computer  solution  is  included  in  the  Computer  Program 
Section  and  provides  for  determining  the  intersecting  value 
of  (4.21a)  and  (4.21b).   It  additionally  calculates  the 
critical  and  amended  loads  for  each  value  of  considered  n. 
These  results  are  summarized  in  Table  I  below. 
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TABLE  I 


n 

V 

v+1 

Pcrit 

P* 

P*/Pu 

0 

0.2500 

1.2500 

l6.936=Pu 

16.936 

1.0 

.1 

0.2368 

1.2368 

15.618 

17.179 

1.014 

.2 

0.2222 

1.2222 

14.288 

17.145 

1.012 

.3 

0.2059 

1.2059 

12.Q44 

16.827 

0.994 

.4 

0.1875 

1.1875 

11.580 

16.212 

0.957 

•  5 

0.1667 

1.1667 

IO.I89 

15.283 

0.902 

.6 

0.1429 

1.1429 

8.758 

14.013 

0.827 

.7 

0.115^ 

1.1154 

7.263 

12.347 

0.729 

.8 

0.0833 

1.0833 

5.656 

10.181 

0.601 

.9 

0.0455 

1.0455 

3.794 

7.210 

0.426 

0.98 

0.0010 

1.0010 

0.507 

1.013 

0.059 

1.0 

0 

— 

0 

0 

0 

The  best  design  exists  between  n=0.1  and  0.2.   Additional 
calculations  indicated  that  n=0.l4  offers  the  largest  buckling 
load  (P*=17.199)  for  the  variable  height  formulation.   This 
design  is  shown  in  Figure  8.   Calculation  of  values  for 
P*/Pu  indicates  that  there  is  little  change  for  the  first  six 
values  of  n.   A  graph  at  P*  versus  n  is  shown  in  Figure  9- 
Figures  27-38  in  Appendix  E  illustrate  the  various  designs 
studied  in  the  interval  of  n=0  to  n=l. 
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Figure  9-   Amended  Load  P*  Versus  Values  of  N,  (case  1). 

C.   EXTENSION  OF  THE  PROBLEM 

The  formulation  for  design  criteria  was  previously  based 
on  the  assumption 


h  =  h  (1  -  2x/i) 


n 


(4.1) 


The  case  will  now  be  considered  for  the  design  based  on  the 
condition 


h  =  h  (1  +  2x/£) 


n 


(4.36) 


This  condition  'is  physically  shown  in  Figure  10.   By  inspec- 
tion of  the  governing  equations  of  the  forementioned  solution, 
it  is  evident  that  the  sign  change  makes  no  difference  in 
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Figure  10.   Variable  Height  Beam  Based  on  h  =  h  (l+2x/£)n. 


the  solution  of  the  problem.   The  identical  critical  loads 
are  again  found.   However,  each  design  is  restricted  to  the 
same  volume  of  material.   This  requires  that  the  midspan 
height  be  decreased  while  the  end  heights  of  the  beam  are 
increased.   This  modification  has  a  significant  effect  on 
the  amended  loads. 

Integrating  (4.36)  and  using  the  same  technique  described 
in  Section  IV-B ,  the  center  height  reduces  to 


h 


n+1 


(2) 


n+1 


-  1 


h 


(4.37) 


while  the  ends  expand  to 


,n+l 


he  "  2(n+l)  lir 


h 


(4.38) 


A  computer  solution  for  this  problem  is  included  at  the  end 
of  the  thesis.   The  results  are  summarized  in  Table  II. 

The  structure  is  less  stable  for  any  comparable  condi- 
tion listed  in  the  previous  section.   Figures  39-48, 
Appendix  E  include  designs  analyzed  for  this  problem.   A 
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TABLE  II 


n 

X 

Pcrit 

P* 

%   Decrease 

0 

1.05851 

16.936 

16.936 

— 

.1 

1.02747 

15.618 

15.023 

11.30 

.2 

0.99220 

14.288 

13.215 

21.97 

.3 

0.95178 

12.944 

11.508 

32.05 

.4 

0.90470 

11.580 

9.891 

41.60 

.5 

0.84907 

10.189 

8.359 

50.64 

.6 

0.78200 

8.758 

6.898 

59.27 

.7 

0.69840 

7.263 

5.490 

67.58 

.8 

0.58920 

5.656 

4.102 

75.78 

.9 

0.43119 

3.794 

2.639 

84.41 

plot  of  the  amended  load  P*  versus  n  is  shown  in  Figure  11. 
Note  that  all  designs  reflect  a  decrease  in  buckling  load 
as  compared  to  Figure  9  and  Table  I.   Also,  there  is  no 
range  of  designs  for  which  the  value  P*  remains  relatively 
stationary . 

D.   LATERAL  BUCKLING  OF  A  CANTILEVER 

A  similar  analysis  [7]  for  a  cantilever  beam  shows  that 
the  buckling- loads  are  determined  by 

J_v  (^7)  =  0  (4.39) 

where  K  is  defined  as 

2 
K  =  P&  (4.40) 

c°c° 

15 
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N 


Figure  11.   Amended  Load  P*  Versus  Values  of  N  (case  2). 

A  computer  solution  was  used  to  calculate  the  first  zero 
of  a  negative  order  Bessel  function  of  the  first  kind,  and 
used  this  result  to  obtain  the  critical  and  amended  loads 
for  the  structure.   These  results  are  summarized  below  in 
Table  III.   The  last  value  is  obtained  by  evaluation  of  the 
Bessel  function  of  the  second  kind  and  use  of  tables  [18] . 
The  values  of  P*  versus  n  are  shown  in  Figure  12.   Note  in 
particular  that  all  values  are  greater  than  the  ultimate 
critical  load,  Pu,  for  the  rectangular  case. 

The  best  design  is  determined  for  the  case  of  n=0.753 
and  is  noted  as  P*=4.9l8.   This  shape  is  shown  in  Figure  13- 
The  range  of  shapes  considered  in  the  formulation  are  shown 
in  Appendix  E. 
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TABLE  III 


n  v  Pcrit         P*         %    Increase 

0  -.2500  4.013=Pu  4.013 

.1  -.2368  3.858  4.239  5.63 

.2  -.2222  3.694  4.433  10.47 

.3  -.2059  3.534  4.595  14.50 

.4  -.1875  3.374  4.724  17.72 

.5  -.1667  3.214  .  4.821  20.13 

.6  -.1429  3.053  4.885  21.73 

.7  -.1154  2.892  4.917  22.53 

.8  -.0833  2.731  4.915  22.48 

.9  -.0455  2.568  4.880  21.60 

1.0    .        0  2.405-  4.810  19.86 


E.   COMPARATIVE  ANALYSIS  OF  SOLUTIONS 

A  comparative  analysis  between  the  solutions  determined 
for  the  cantilever  and  simply  supported  beams  provides  some 
interesting  insight  into  the  actual  physics  of  the  problem. 
The  results  obtained  from  the  cantilever  solution  appear  to 
be  a  near  optimum  shape  for  the  structure.   At  first  sight 
analysis  of  the  simply  supported  beam  does  not  appear  to 
yield  the  optimum  shape,  although  a  better  design  than 
that  of  classical  solution  is  obtained.   However,  a  close 
examination  of  the  physics  of  the  problem  leads  to  a  prediC' 
tion  for  the  possible  optimum  shape,  although  it  was  not 
analytically  determined. 
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Figure    12 


Amended  Load  P*  Versus  Values  of  N 
for  Cantilever  Beam. 


The  value  n=0.l4  provided  a  best  buckling  load  of 
P*  =  17.199  for  the  simply  supported  beam.   This  corresponds 
to  a  gain  of  only  1.55*  over  that  of  a  constant  rectangular 
cross  section.   For  the  cantilever,  it  was  found  that 
P*  =  4.918  (n=0.75)  which  is  a  22.5??  gain  over  that  of  the 
rectangular  section.   Physically  these  results  can  be 
explained  when  the  lateral  bending  and  twisting  moments  of 
the  two  problems  are  examined. 

Recalling  equations  (3.4)  and  (3.6)  for  the  rectangular 
beam  yields  the  following  principle  axes  moments 


P  ,  \    j.    p  (l 

ME  =  ~  2  ^ul~u-)    2  (2~ 


x)  u 


(4.41) 
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Mn  =  |  (f  -  xH  (4.42) 

Evaluation  at  the  boundary  condition,  x=0 :   u  =u  and  u'=0 
shows  that  the  twisting  moment  Mr=0  and  bending  moment  M 


assumes  its  maximum  value.   At  the  ends  of  the  structure 


x  =  p-;  <J>  =  0  and  u=0.   This  implies  that  bending  moment  M  =0 
and  twisting  moment  M,.  is  a  maximum. 

The  combination  of  these  two  moments  appear  to  restrict 
the  beam  from  assuming  an  optimum  shape  which  differs 
significantly  from  the  uniform  case.   In  fact,  if  the 
structures  depicted  in  Appendix  E  are  examined,  it  may  be 
noted  that  as  the  beam  approaches  a  pyramidal  shape,  the 
less  stable  the  structure  becomes.   Physically,  this  seems 
to  imply  that  the  vanishing  end  area  offers  less  and  less 
resistance  against  the  twisting  moment,  and  the  structure 
will  readily  buckle. 

It  is  apparent  and  expected  that  the  existing  moment 
conditions  govern  the  optimum  shape  of  the  structure. 
Physical  intuition  would  indicate  that  the  maximum  amount 
of  material  be  located  where  the  maximum  moment  existed. 
If  M   and  Mr  were  equal  and  a  linear  variation  could  be 
assumed,  the  resulting  optimum  shape  would  actually  be 
that  of  a  rectangular  beam.   For  the  problem  considered, 
M   appears  to  have  a  domineering  influence  on  the  optimum 
design.   If  this  was  not  the  case,  then  the  solutions 
determined  in  Section  IV-C  would  have  indicated  a  "best" 
buckling  load  greater  than  those  for  the  present  solution. 
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The  actual  moment  values  for  M   and  Mr  are  not  known 
and  cannot  be  determined  unless  the  displacement  and  twist 
functions  u  and  <J>  are  known  exactly.   Additionally,  these 
functions  are  coupled  to  the  variable  cross  sectional  area 
as  indicated  by  equations  (3-3  a,b,c).   Analyzing  equations 
(4.41)  and  (4.42)  and  Figure  4  leads  to  the  following  inter- 
pretation.  The  moment  representation  of  M..  appears  to  be 
parabolic  while  that  of  M   (based  on  a  parabolic  <J>)  appears 
to  be  cubic  and  exhibits  an  inflection  point.   These  patterns 

are  shown  in  Figure  14.   The  maximum  combined  moment  must 

o 

then  occur  between  x=0  and  x=  p-  . 


:•: 


M   =  lateral  bending  moment 
M^  =  twisting  moment 


PREDICTED  OPTIMUM  SHAPE 


♦-  x 


x=0 


x=£/2 


Figure  14. 


Expected  Moment  Configuration  and 
Approximate  Optimum  Design. 


On  the  basis  of  moment  evaluation  alone,  the  optimum 
design  would  apparently  be  of  a  form  similar  to  the  shape 
shown  in  Figure  14.   Since  M   is  expected  to  be  greater 

than  M,_  the  design  would  lump  the  largest  amount  of  material 

I 
somewhere  between  x=0  and  x=  p-  but  closest  to  x=0. 
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Available  time  has  prevented  analysis  of  this  expected 
"optimum"  shape.   A  Rayleigh  Ritz  formulation  for  this 
design  was  attempted  by  assuming  the  displacement  and  twist 
functions  and  using  a  Lagrange  interpolating  polynomial  to 
approximate  the  design.   The  results  indicated  a  buckling 
load  slightly  less  than  that  of  the  uniform  rectangular 
case . 

Because  of  the  low  gain  in  efficiency  (1.55$)  from  the 
best  design  obtained  from  the  variable  height  solution,  it 
appears  that  the  true  optimum  shape  will  not  reflect  a 
large  performance  gain  over  that  of  the  rectangular  case. 
In  fact,  the  determined  value  from  the  variable  height 
solution  is  expected  to  be  very  close  to  the  actual  optimum 
design. 

F.   VALUES  OF  N  GREATER  THAN  ONE 

The  analyses  of  the  problems  discussed  have  been  restric- 
ted to  the  range  for  n=l.   If  n  is  greater  than  one  the 
actual  curvature  of  the  beam  design  becomes  concave.   For 
the  simply  supported  beam  this  physically  Implies  that  the 
structure  must  become  less  stable  than  the  limiting  pyramidal 
shape.   This  is  evident  since  even  less  area  is  distributed 
toward  the  beam  ends.   But,  it  was  found  earlier  that  the 
pyramidal  shape  is  actually  unstable  for  any  load.   By  this 
reasoning  the  problem  must  be  restricted  to  the  range  of 
values  used  to  have  physical  meaning. 

Mathematically,  if  the  value  of  n  is  allowed  to  be 
greater  than  one,  then  equation  (4.15)  changes  to  the  form 
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1-n 
4.  =  Y  2   [a1J_v(p)  +  ct2Jv(p)  ]  (4.45) 

for  which  a,  must  be  zero.   The  resulting  roots  introduce 
a  singularity  condition  and  the  solutions  are  misrepresen- 
tative  of  the  problem.   Consequently,  it  is  deemed  necessary 
for  the  value  of  n  to  be  restricted  to  the  interval 
0  <  n  <  1  to  provide  meaningful  results. 
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V.   CONCLUSIONS 

For  the  simply  supported  beam  the  assumed  power  law 
restricts  the  number  of  designs  which  may  be  analyzed. 
However,  a  good  indication  of  the  problem  behavior  can  be 
determined  from  the  results.   The  number  of  designs  which 
remain  relatively  stable  do  not  differ  radically  from  the 
uniform  case.   As  the  designs  converge  to  a  pyramidal 
shape  they  rapidly  become  less  stable. 

It  appears  that  the  combination  of  the  two  moments , 
M   and  M*.  govern  the  design.   As  the  beam  approaches  a 
pyramidal  shape,  less  area  becomes  available  to  resist  the 
twisting  moment  and  the  beam  readily  buckles.   Thus  the 
interaction  of  the  two  moments  restrict  the  beam  from 
attaining  a  shape  which  differs  significantly  from  the 
rectangular  case. 

The  effect  of  the  moment  values  on  beam  design  is 
better  demonstrated  from  the  cantilever  results.   At  the 
free  end  of  the  beam  there  is  no  moment  and  the  design  calls 
for  very  little  material  to  be  located  at  that  position. 
The  interaction  of  the  two  moments  seem  to  produce  a  maximum 
slightly  away  from  the  fixed  end  -  otherwise,  the  design 
for  the  best  buckling  load  would  have  indicated  the  pyra- 
midal shape.   The  cantilever  solution,  thus,  helps  to 
support  the  validity  of  the  simply  supported  beam  results . 
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The  true  optimum  shape  for  the  simply  supported  beam 
has  not  been  determined  by  the  variable  height  formulation 
However,  the  expected  optimum  shape  was  approximated  from 
analyzing  the  physics  of  the  problem.   This  shape  is  not 
expected  to  offer  any  significant  increase  over  the  design 
obtained  from  the  variable  height  formulation. 
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APPENDIX  A 
A  STRUCTURAL  OPTIMIZATION  PROBLEM 

An  example  of  a  structural  optimization  problem  which 
illustrates  the  usefulness  of  the  variational  method  of 
Section  II  is  that  of  obtaining  the  optimum  shape  of  a  uni- 
formly loaded  beam.   This  problem  is  a  "strength"  problem 
in  contrast  to  the  buckling  problem  which  exhibits  no 
strength.   For  the  problem  considered,  bending  is  in  the 
vertical  plane  of  symmetry  (Fig.  15).   The  height  of  the 


Figure  15-   Variable  Width  Uniformly  Loaded  Beam. 

beam  is  taken  as  constant,  and  a  uniform  load  is  applied. 
1 .   Statement  of  the  Problem 

The  design  b(x)  is  to  be  determined  for  a  variable 
width  which  provides  a  maximum  strength  capacity  for  the 
structure.   Unless  some  specific  constraints  are  imposed  on 
the  problem,  results  may  well  be  meaningless.   For  instance, 
if  a  minimum  volume  structure  was  considered  with  no  other 
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conditions  than  equilibrium  and  boundary  conditions  to  be 
satisfied,  then  it  is  entirely  reasonable  to  deduce  a  zero 
volume  structure  with  infinite  strength  capacity  [6]  . 

For  the  strength  problem  the  primary  constraint  is  that 
the  material  must  not  fail.   Thus,  for  an  elastic  structure 
composed  of  material  of  finite  strength  capacity,  the 
specific  strain  energy  (ij;)  must  be  no  greater  than  a  speci- 
fied maximum.   This  maximum  is  denoted  as  ^,  such  that, 

*  i  £  (A.l) 

at  any  point  in  the  structure.   The  inequality  of  (A.l)  may 
be  written  as  an  equality  constraint  by  introducing  a  scale  r 
variable  y(x)  to  take  up  any  existing  slack  in  the  function. 
Using  the  "slack  variable"  gives 

i>   +  y2  =  i>   .  (A. 2) 

For  the  simply  supported  beam  the  moment  of  inertia  with 

1       3 
respect  to  the  principle  bending  axis  is  I  =  -^p   b(x)hn   or 

I  =   CA(x) .   The  strain  energy  per  unit  length  at  any  posi- 

2 
tion  along  the  beam  is  known  to  be  l/2ECAv"   where  v  repre- 
sents the  displacement  function.   If  the  strain  energy  per 
unit  length  is  divided  by  the  cross  sectional  area,  A,  at 
any  position,  then  specific  strain  energy  per  unit  area  per 
unit  length  may  be  written 

if;  =  l/2ECv"2.  (A. 3) 

An  energy  constraint  equation  for  the  problem  is  then 
determined, 
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J  -  (|  ECv"2  +  y2)  =  0  .  (A. 4) 

A  second  constraint  on  the  problem  restricts  the  beam  to  be 

designed  from  a  given  volume  of  material  V  .   This  is  the 

isoperimetric  constraint  introduced  in  Section  II.   Thus 

I 
f      A(x)dx  =  V   .  (A. 5) 

0  ° 

The  total  potential  energy  of  the  system  may  be  written 

T  ■  /[|  ECAv"2  -  Pqv]  dx,  (A. 6) 

where  the  first  term  represents  the  strain  energy  due  to 
bending  and  the  second  term  defines  the  external  work  of 
the  system  due  to  p  [10] . 

The  volume  constraint  (A.  5)  and  energy  constraint  (A. 4) 
are  now  accommodated  through  the  formation  of  the  augmented 
functional  T* , 
I 
0 


T*  =  /   [fecAv"2-  p  v  -  A,A-A0(>  -i-ECv,,2+Y2)]dx  (A. 7) 

n     &  O      -L     c      d 


where  A,  and  A„  are  Lagrangian  multipliers. 

Performing  a  variation  with  respect  to  v,  A,  and  y 
yields  three  Euler  equations. 

6  T*=0:   [ECAv"  +  ECv"A„]  -  p   =0  (A. 8) 

v        L  2J    Ko 

6AT*  =  0:   3gECv"2  -  X±   «  0  (A. 9) 

6  T*=0:   2A2y  =  0  .  _         (A. 10) 

The  natural  boundary  conditions  are  v(0)=v(L)=0  and 
M(0)=M(L)=0. 


58 


Integrating  equation  (A. 8)  twice  yields 

EC(Av"  +  A2v")  =  h   pQX2  +  C-jX  +  C2  .  (A. 11) 

"po 
for  which  the  above  boundary  conditions  determine  C,  =  — r— 


and 


Cp  =  0.   Consequently  (A. 11)  becomes 


v"(A+X2)  =  2g§  (X  -£x).  (A. 12) 

Using  equation  (A. 9)  gives 

2  A  h 

v"  =  (— i)  (A. 13) 

EC 

which  physically  implies  that  the  deflection  curvature  is 
constant . 

From  equation  (A. 10)  there  exist  three  possibilities 

i)   X2  E  0,y  E.O  (A. 14) 

ii)   A2  %    0,y  =  0  (A. 15) 

iii)   A2  e  0,y  ^  0  .  (A. 16) 

Assuming  y  =  0  and  examining  equation  (A. 4),  it  may  be 
observed  that 

v"  =  (§§)  (A. 17) 

Comparing  this  with  (A. 13)  implies 

J  =  X.  .  (A. 18) 

Recalling  [11]  for  a  simply  supported  beam  under  a 
uniform  load  Po,  the  moment  at  any  position  X  along  the  beam 
may  be  represented  by 

M(x)  ■  -y  (X^  -  £x)  (A. 19) 
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then  equation  (A. 12)  may  be  rewritten 

v"(A+X2)  =  j^j  M(x)  .  (A. 20) 

But,  the  moment  is  also  known  as  M=ECAv" ;  consequently, 

v"(A+X2)  =  v"A  .  (A. 21) 

From  the  above  equation  it  is  obvious  that  X~  =  0  and  (A.l4) 
is  satisfied.   This  condition  implies  that  the  energy  formu- 
lation actually  rejects  an  energy  constraint.   This  can  be 
observed  by  examining  equation  (A. 8),  which  should  be  the 
equilibrium  equation  of  the  system.   The  condition  is  valid 
only  if  X_  is,  in  fact,  zero.   The  rejection  of  an  energy 
constraint  appears  to  be  inherent  to  the  energy  formulation. 

This  formulation  seeks  to  maximize  the  maximum  stress 
at  every  position  x  (i.e.  strain  energy  density)  which  is 
different  than  maximizing  specific  strain  energy  as  was  for- 
mulated in  the  problem. 

Therefore,  equation  (A. 21)  becomes 

P     n 
v"A(x)  =  £e§  (x  -Jtx)  (A. 22) 

from  which  the  variable  area  may  be  determined 

A(x)  =  P°(x  ~£x)  (A. 23) 


(8  $  EC)51 


Substituting  for  A(x)  into  the  isoperimetric  equation  (A. 5) 
and  integrating,  the  load  Po  is  found 

po  =  _  6Vo  (8^EC)a-2  (A.24) 

I5 
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The  result  determines  the  optimum  shape  of  the  cross  sec- 
tional area,  or  in  terms  of  the  width 


b(x)  = 


6Vo 
£3h 


(ix-xd)   . 


(A. 25) 


a  = 


Since  the  stress  at  any  section  may  be  represented  as 

Mc 
=  -y-,  substitution  for  M  and  c  (c  =  h  /2)  yields 


a  = 


Eh    oT  h 
2   ^EC; 


Eh    „,  3* 

o  /2Xv  2 

T"  (EC}   =  a 


(A. 26) 


This  represents  the  maximum  allowable  stress  at  any  section 
along  the  beam.   It  is  noted  additionally,  that  the  stress 
is  independent  of  any  position  x  along  the  beam.   Thus,  the 
stress  is  maximized  at  the  outer  fibers  (h  /2)  for  any 
position  with  respect  to  the  optimum  shape.   The  optimum 
shape  is  noted  to  be  parabolic  and  is  shown  below  (Pig.l6) 
with  stresses  indicated.   The  design  is  comparable  to  the 
uniform  strength  design  listed  by  Esbach  [20] . 


o^a- 


bOO-f^Ux-x4) 


Fierure  16.   Optimum  Shape  for  Variable  Width  Beam  Design 
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APPENDIX  B 
EFFECT  OF  THE  LOAD  P  AT  A  DISTANCE  FROM  THE  CENTROID 

In  Section  III  the  assumption  was  made  that  the  load  P 
should  be  located  at  the  centroid  of  the  beam.  The  effect 
of  the  geometrical  location  of  the  load  is  significant  in 
the  energy  formulation  analysis.  For  the  lateral  buckling 
case,  the  external  work  associated  with  the  unstable  state 
is  represented  by 

1/2 
p  /    4>uM  (1/2    -  x)  dx  (3.20) 

0 

This  condition  results  from  analyzing  a  longitudinal  element 
dx  at  point  D  in  Figure  4.   Considering  bending  of  this 
element  in  the  £C  plane  with  the  cross  section  mn  assumed 
fixed,  the  end  of  the  beam  then  describes  an  infinitely 
small  arc 

u"U/2  -  x)  dx  (B.l) 

in  the  ££  plane,  for  which  the  vertical  component  may  be 
represented 

<$>u"(l/2    -  x)  dx  (B.2) 

The  external  work  is  P  times  the  integral  of  this  term  over 
the  area  considered. 

Now  consider  the  case  where  the  load  P  is  located  at  the 
top  of  the  beam  (Figure  17) •   The  external  work  must  now  be 
determined  from  lowering  of  the  load  by  amounts  A,  and  A2 

[191. 
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h/2 


Figure  17.   Effect  of  the  Load  P  Located  at  Top  of  Beam. 


A,  is  simply  lowering  of  the  load  due  to  the  change  in 

centroidal  height  during  lateral  buckling  and  is  the  same 

as  (3.20).   The  component  A„  represents  the  lowering  due  to 

twisting  of  the  beam  through  an  angle  <t>  and  can  be  repre- 

p 

sented  as   h/2  (1  -  Cos  <$>)    -   ^~  .   Thus  an  additional  term 


mu 


st  be  considered  if  the  load  is  placed  away  from  the  beam 


Ph   2 
centroid.   This  term  reduces  to  -jp  4>  . 
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APPENDIX  C 
TRANSFORMATION  OF  EQUATION  (4.11) 


dfi+  n   di+  R2  2(1-*)    =  Q  (1|>n) 


dy2   Y   dY 


P  =  o Y 

2-n 


K_  v2-n  (4.12a) 

-n 

1-n 
4>  =  Y  2   6(Y)  (4.12b) 

from  (4.12b) 

-1-n      1-n 

H7  "  ~Y    6  Y   H7  (c<1) 

where 


d$    dB_  dp_ 
dy    dp   dy 

Substituting  to  get 


(C.2) 


-1-n      3/  -,  _  \ 

di  _  lzH  v   2       2U  nj    d§_  r         , 

H7  "  —  Y     6+Y       K^-  (C.3) 

Differentiating  again 

p  -  3-n  -l-3n 

O.  _  (1-n)   (-1-n)     2      pf    .     2   d§_ 
dy2 ~    2     Y      6  +  2(l-n)K  Y      ^ 

.  ^2  5/2  (1-n)  d26  (C.4) 

+  is.  y  — p 

Substituting  from  (4.12b),  (C.3)  and  (C.4)  into  (4.11)  and 
arranging  in  tabular  form: 
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dfi 
dpi 

K2  !(l-n) 

K  Y 


dB 
dp 


l-3n 


2(l-n)KY 


l-3n 


nKY 


-3-n 
n(n-l)     2 


v2  !(l-n) 

K  Y 


d2(j) 
dY2 


n  d$_ 
Y  dY 


K2  2(l-n) 


|(l-n)  ,2 
K  Y 


(l-3n) 


M-  +  (2-n)KY   2 


dp' 


dp 


n   +  n    1W 


T~ 


IT- 


d      +  K  Y 


8     (C5) 


The  operator  frorr  equation  (C.4)  is  shown  at  the  right  side 

of  the  table. 

Substituting    (il.l2a)    and   letting 


1-n 
v   =    2T2^HT 


(C.6) 


Equation  (5)  becomes 


S-i+  i  ||+  (1 

dp2     p   dp 


%)  6 


=  0 


ft. 13) 
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APPENDIX  D 
BESSEL  EQUATION  CURVES  FOR  LATERAL  RUCKLING 


-JKL 


CASE:  N*0 


rf-SCRLF.----l-.GOE:.+00  UNITS  \KZH. 
^-SCRLF.---2.00EhOG  UH]TS  ]NCH- 

R.    Ln    BURNS   LCDR,    U5N 

BESSEL   EQ-    CURUES   FOR   BUCKLING 


Figure  18. 
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I  '  .r 


IrJKl 


K    '^i-.: 


n  n  (7 


CASS:  N=0.1 


F.+  OG   UNITS    INCH, 
^-SCRLF--?.aOE+OG  UNITS    INCH. 

R,    l-    BURNS   LCDR,    USH 
BESSEL    EQ-    CURUE5   FOR   BUCKLING 


Figure  19«. 
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FJK 


JW 


CASE:    11=0.2 


K-SCRLF----J..OQE+00   UNITS    INCH- 

*-SCRLE-2.00E-*0a   UNITS    INCH. 

R3    L-    BURNS   LCOR,    USN 
BESSEL   EQ-    CURUES   FOR   BUCKLING 


Figure  20. 
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K-5LRLF.-1.0DE+D0  UHJTS    INCH,  case:  k=o.3 

V-5CRLF-2. QGE+OO  UNJT5    INCH. 

R.    L.    BURNS   LGDR,    USN 
BESSEL   EQ.    CURUE5   FOR   BUCKLING 

Figure  21. 
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K-SCRLF-1.00E-+00  UNITS  INCH. 
?-5CRLE«2.00E+00  UNITS  INCH. 


CASS:  N=0.4 


R.    L.    BURNS   LCOR,    USN 
BES5EL   EQ.    CURUES   FOR   BUCKLING 


Figure  22. 
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K-5tflLE^1.00E+0G  UHJT5  INCH. 
?-5CnLE.---2.  OCE+OG  UNJT5  INCH. 


CASS:  N=0.5 


R.    L.    BURNS   LCDR,    USN 
QES5EL    EQ.    CURUE5   FOR   BUCKLING 


Figure  23. 
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K-SCRLE^-l.'OOE-vOa   UNITS    INCH.  case:  k=o.6 

H,-SCRLF-2.00E-+DQ   UNITS    INCH. 

R.    L.    BURNS   LCQR,    USN 
BESSEL   EQ.    CURUE5   FOR   BUCKLING 

Figure  2k. 
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rj<n 


K-5CnLr-l.-0Dr.-0D   UNITS    INCH.  case:  n=o.7 

^j-5CRlP--5.D0L-DD   UHJTS    INCH. 

"  R,  [_■  BURNS  LGDR,  U5N 
BESSEL  EQ-  C'URUES  FOR  BUCKLING 

Figure  25-. 
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^-5CRLE---5.0CE-tOO  UNITS   INCH. 

R.    L.    BURNS   LCQR,    USN 
BE5SEL   EQ.    CURUES   FOR   BUCKLING 

Figure  26. 
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Figure  2  7., 
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APPENDIX  E 
VARIABLE  HEIGHT  BEAM  DESIGNS 
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COMPUTER  PROGRAMS 
1.   BESSEL  EQUATION  CURVE  GENERATION 
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4.       CANTILEVER   BEAM 
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5.       GENERATE    BEAM   DESIGNS 
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